Casimir— Polder forces on excited atoms in the strong atom— field coupling regime 
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Based on macroscopic quantum electrodynamics in linear media, we develop a general theory of 
the resonant Casimir-Polder force on an excited two-level atom in the presence of arbitrary linear 
bodies, with special emphasis on the strong-coupling regime where reabsorption of an emitted photon 
can give rise to (vacuum) Rabi oscillations. We first derive a simple time-independent expression for 
the force by using a dressed-state approximation. For initially single-quantum excited atom-field 
systems we then study the dynamics of the force by starting from the Lorentz force and evaluating 
its average as a function of time. For strong atom-field coupling, we find that the force may undergo 
damped Rabi oscillations. The damping is due to the decay of both the atomic excitation and the 
field excitation, and both amplitude and mean value of the oscillations depend on the chosen initial 
state. 

PACS numbers: 12.20.-m, 42.50.Vk, 42.50.Nn, 32.70.Jz 



I. INTRODUCTION 

It is well known that the presence of macroscopic bod- 
ies can drastically change the properties of the electro- 
magnetic field compared to those observed in free space. 
A prominent manifestation of the body-induced change of 
the ground-state fluctuations of the field is the Casimir- 
Polder (CP) force experienced by an atom when placed 
within an arrangement of unpolarized ground-state bod- 
ies (for a review, see Ref. [l|). Originally, CP forces were 
studied for ground-state atoms where the CP potential 
can be identified with the position-dependent part of the 
atom-field coupling energy Q • Since in this case the cou- 
pling involves only off-resonant, virtual transitions of the 
system, the coupling energy can be calculated by means 
of time-independent leading-order perturbation theory. 

For excited systems, on the contrary, real transitions 
must be taken into account. In particular, when applying 
the method to an atom in an excited energy eigenstate 
which interacts with the body-assisted electromagnetic 
vacuum, one finds that the corresponding CP potential 
can be significantly enhanced The enhancement is 
due to the influence of now possible real transitions to 
lower states, with the dependence of the force on the 
atomic transition frequencies exhibiting the typical dis- 
persion profiles in the vicinity of medium resonances. 
Since these transitions are also responsible for the de- 
cay of the atomic excitation, the application of the static 
approach to the CP force on an excited atom becomes 
questionable. 

Moreover, perturbative methods are known to break 
down when an atom near-resonantly interacts with a 
body-assisted narrow-band field such that the strong- 
coupling regime is realized; this is typically the case when 
the bodies form a resonatorlike structure. For a two-level 
atom that (according to the Jaynes-Cummings model 
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[H) is assumed to near-resonantly interact with a single 
monochromatic mode, it is well known that the energies 
of the exact excited energy eigenstates of the coupled 
system — the so-called dressed states — are symmetrically 
shifted above/below the mean of the unperturbed energy 
eigenvalues. As the dressed-state energies depend on the 
atomic position, they can be interpreted as CP potentials 
0, Q — in generalization of the static approach based on 
leading-order perturbation theory. In particular in the 
case of a perfect standing wave in a cavity, it turns out 
that, depending on the upper/lower dressed state the sys- 
tem is prepared in, the atom is repelled from/attracted 
to the antinodes of the wave — an effect that can strongly 
affect the motion of the atom 0] • Calculations imply that 
atoms sufficiently slowly incident on a cavity in a direc- 
tion parallel to the standing wave can be reflected from 
the cavity due to the influence of the upper dressed state 
(fill, Q. Here, the dependence of the transmission prob- 
ability on the velocity of the incident atoms displays a 
resonance structure [lfj which may be used to construct 
a velocity filter for atoms [ll|> [12] • In the case of atoms 
incident in a direction normal to the standing wave, the 
atoms may be deflected during their passage through the 
cavity 

In fact, even if a dressed state was prepared very accu- 
rately, the system would not stay in this state forever, but 
it would undergo a temporal evolution due to unavoid- 
able dissipation, and so would the CP force. To describe 
the time-dependence of CP forces acting on excited sys- 
tems, a dynamical approach is required. For weak atom- 
field coupling, it has been shown that when an initially 
excited atom exponentially decays to the ground state, 
the associated CP force shows a similar transient behav- 
ior and eventually changes to the ground-state force in 
the long-time limit p| [l5[ • A two- level atom which is 
initially prepared in the upper state and which strongly 
interacts with a (vacuum) cavity field undergoes Rabi os- 
cillations which are damped due to dissipation; this effect 
may be phenomenologically incorporated in the Jaynes- 
Cummings model via damping constants fl6j . As simu- 
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lated in Refs. [IT], [l8|, the losses can be compensated for 
in a controlled way by introducing an external pump- 
ing laser. Furthermore, by appropriate choice of the laser 
intensity and frequency, it can be ensured that the sys- 
tem remains in a certain steady state associated with an 
effective CP potential, both the sign and magnitude of 
which can be continuously controlled. Such a setup has 
been used to trap single (cold ) at oms in the antinodes 
of a standing wave in a cavity [H, US HH, HH ■ Since the 
coupling of the atom-cavity system to the pump laser de- 
pends on the atomic position within the standing wave, 
it is possible to continuously monitor the atomic motion 
by recording the intensity of the pump laser leaving the 
cavity [IS El, US Hi, 25]. This information can be used in 
a feedback mechanism to enhance the trapping efficiency; 
whenever the atom tends to leave the wave antinode, the 
trap de pth is enlarged by increasing the pump laser in- 
tensity in. 

In the present paper, we study the resonant CP force 
exerted on a two-level atom that strongly interacts with a 
body-assisted electromagnetic field in more detail, where 
the initial state may be a (coherent) superposition of (i) 
the atom being in the upper state and the body-field 
system being in the ground state and (ii) the atom be- 
ing in the lower state and the body-field system being in 
a single-quantum excited state. On the basis of macro- 
scopic quantum electrodynamics (QED) in linearly re- 
sponding media, we aim at generalizing the theory with 
two respects. Firstly, we allow for an arrangement of bod- 
ies of arbitrary shape and material which are character- 
ized by their space- and frequency-dependent complex 
permittivity and permeability, thus accounting for both 
material absorption and dispersion in a natural way. The 
spectral and spatial structure of the body-assisted field 
follows from the Green tensor of the associated macro- 
scopic Maxwell equations, thus generalizing the idealized 
case of a single perfect standing wave in a cavity to more 
realistic systems. Secondly, we include the nontrivial tem- 
poral evolution of both the state of the system and the 
CP force in the theory by developing a dynamical ap- 
proach to the CP force. 

The article is organized as follows. After an introduc- 
tion of the system Hamiltonian within the framework of 
macroscopic QED (Sec. [TTJ) , we first develop a static ap- 
proach to the CP force by assuming a dressed-state-like 
preparation (Sec. IIII A[) . We then develop a dynamical 
theory of the CP force (Sec. IIII B|) by considering the tem- 
poral evolution of the quantum averaged Lorentz force. 
We end with a summary in Sec. IIVI 



II. BASIC EQUATIONS 

A. Sketch of the quantization scheme 

Consider a neutral atom consisting of particles of 
charges q a , masses m a , positions f a , and canonically con- 
jugate momenta p Q (a = 1,2, . . .) which interact with the 



electromagnetic field in the presence of linearly respond- 
ing bodies. In particular, let us consider an arrangement 
of magneto-electric bodies, which is described by the 
spatially varying, complex permittivity e(r, uj) and per- 
meability ^t(r, w). In electric-dipole approximation, the 



multipolar-coupling Hamiltonian that governs the dy- 
namics of the system can be given in the form [l|, IbH . fl5| 

H = H A + H F + H AF . (1) 



Here. 



is the atomic Hamiltonian, with 



Pa(i-) = y^q a i a / dA<5(r - f A - Ar Q ) (3) 



being the atomic polarization relative to the center of 
mass 



r Q 



(dia = J2 a m a)i where 

r a = r a — r A 



(4) 



(5) 



are the relative particle coordinates. The Hamiltonian 
H-p of the combined system of the electromagnetic field 
and the magneto-electric bodies, and the atom-field 
(dipole-)coupling term H A f can be expressed in terms 
of Bosonic (collective) variables f\(r,u>) and f J(r, oj), 

[h(r,u;)Jt l (r',oj')]=S X yS(u;-u;')S(r-r'), (6) 
[f A (r,w),f v (rV)] =0 (7) 

(A, A' G {e, m}), as follows: 

/poo 
d 3 r I d^fj(r,^)-f A (r,^), (8) 

. \ — i .in ^ 



H AF = -d-E(f A ), 



(9) 



where 



a. a. 

is the electric dipole moment of the atom, and 



E(r)= J dLJ y"dVG A (r,r>)-f A (r>) + H. 

A— e,m 



(ii) 
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with G\(r, r',ui) being related to the classical (retarded) 
Green tensor G(r, r', u>) according to 



G e (r, r', w) = i J Im e(r>, u>) G(r, r',u>), (12) 



CT V 7T£o 



G m (r, r',u)=i- J Im n(r', lu) [ V'x G(r', r, w)] T . 

C V 7T£0 

(13) 

Note that E(r) has the physical meaning of a displace- 
ment field with respect to the atomic polarization, and 
the body-assisted induction field reads 



«w = E / *v /. S 



V x «(r, V x e(r, u>) 



G{y,y',uj) = <5(r-r') 



(15) 

[«(r, uj) = fj. 1 (r,cj)] together with the boundary condi- 
tion 



G(r, r', u>) — > for |r — r'| — > oo. 
It has the following useful properties [2^ |: 
G*(r,r» = G(r,r',-u;*), 

G(r,r',w) = G T (r',r,w), 



(17) 
(18) 



J < ^ s {~2 Im e ( s : w ) G ( r > s > w ) ' G*(s, r', w) 

- Im «(s, w) [V s x G(s, r, w)] T [V s x G*(s, r', w)] } 

= ImG(r,r',w). (19) 
Combination of Eq. HU) with Eqs. |12|) and ((13J) leads to 

^ / 'd 3 S Gx(r,s, W )-Gl T (r',s, w ) 

A— e,m 

= ^w 1 ImG(r l r , 1 w), (20) 

For an extension of the quantization scheme to arbitrary 
linear bodies (including nonlocally responding ones), we 
refer the reader to Ref. 1281. 



B. Atom— field coupling 

Let us assume that a single atomic transition is cou- 
pled near-resonantly to a body-assisted electromagnetic 



narrow-band field. In this case, it is appropriate to use 
the model of a two-level atom, where only two atomic 
energy eigenstates, denoted by |1) and |0), are involved 
in the atom-field interaction, so that the atomic Hamil- 
tonian ([2]) effectively reduces to 



H A = hhu w &, + ±(£ + Ei)I 



(21) 



[lj w = (E 1 -E )/H; <r z = |l)(l|-|0)(0|; I: unit operator], 
and the electric dipole moment (fTQ|) takes the form 



(22) 



d = d i<7 + H.c. 



(d mn = (m|d|n), d mm = 0, <7=|0)(1|). Consequently, upon 
recalling Eq. (TTT]) , the interaction Hamiltonian © can be 
written as 



x VxG A (r,r»-f A (r»+H.c. (14) jj 

For magneto-electric bodies, the Green tensor is de- 
fined by the differential equation 



AF 



^ J dcvjd 3 r (dioir 1 + d i<r) -G A (rA, r, w)-f A (r, w) 

\— c,m ^ 

+ H.c. (23) 



For the following, it is useful to introduce position- 
dependent photon-like creation and annihilation opera- 
tors a(r,ui) and a'(r,ui) according to the definition 



(16) a(r,w) 



hg(r 



L^V / dVd 10 .G A (r,r».f A (r», 



where 



K r , w ) = -\/^^ 2 dio-ImG(r,r,w)-d 01 



(24) 



(25) 



Substituting Eq. (|24|) into Eq. ([23|) and restricting our 
attention to the rotating-wave approximation, we may 
write the interaction Hamiltonian in the form 

/>oo 

H AF =h dujg(r Al u)a(r A , uj)a ] + H.c, (26) 
Jo 

showing that <7(r, u>) may be regarded as a generalized 
atom-field coupling strength. 

The commutation relations of a(r,u>) and a'(r,u;) can 
be found from the commutation relations © and Q by 
making use of the definitions (|2"4")l and (|23|) and the inte- 
gral relation (12TJ)) , resulting in 

[a(r )W ),ot(r',a/)] = J^f) ^ {w-u/), (27) 



g(r,uj)g(r',cj) 



[a(r,u}),a(T,w')] = 0, 
with 



(28) 



g(r,r',uj) = £ w 2 d 10 -ImG(r,r',u;)-doi. (29) 
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Note that 

[a(r,u),tf(r,w')] =S(uj-u'). (30) 

The definition of the ground state |{0}) of the system 
composed of the bodies and the electromagnetic field, 
f A (r,w)|{0}) =0 (VA,r,w), implies that 

a(T,w)\{0}) = Vr,w. (31) 

The operators a'(r, u) can be used to define single- 
quantum excitation states 

|r,o,)=a+(r,a,)|{0}>. (32) 

From Eq. (|27f it then follows that these states are or- 
thogonal with respect to frequency, but not with respect 
to position, 

(r, u\r', u') = fi^TT- \ S ^ ~ ^ ^ 
g(r,uj)g(r',oj) 

for equal positions we thus have 

(T,w\r,w') = 5{w-u;') (34) 

[cf. Eq. (|30"|) ]. Furthermore, the states \r,ui) are eigen- 
states of Hp carrying one quantum of energy Hu>, 

H F \r,u) = tkj\r,Lo), (35) 

as can be seen from the commutation relation 

[H F) a}(r,u;j\ = tbuatfaw), (36) 

which is a consequence of Eqs. j8|) and (f24|) together with 
the commutation relations ^ and ([7]). 

Equation ([26]) implies that the states |r A ,u;) (w~u 10 ) 
as given according to Eq. (|52"|) describe the single- 
quantum excited states of the body-assisted electromag- 
netic field which may be obtained after an atom posi- 
tioned at fa has undergone an (electric-dipole) transition 
|1) — > |0). Hence, these states can be regarded as span- 
ning the state space of all the states that are resonantly 
coupled to an initially excited two-level atom via Haf as 
given by Eq. (|26[) . which is why it will be sufficient to 
work with them in the following. 



III. THE CASIMIR-POLDER FORCE 

As shown by Casimir and Polder Q, dispersion forces 
on ground-state systems can be obtained by means of 
time-independent perturbation theory, by calculating the 
shift of the respective (unperturbed) ground-state energy 
due to the atom-field coupling for given atomic center- 
of-mass position, 

AE = A£(0) + u(r A ), (37) 



and identifying its position-dependent part U(r A ) as the 
potential from which the force can be derived, 

F(r A ) = -V A U(r A ). (38) 

When dealing with a system prepared in an excited state, 
this simple approach may fail for two reasons. Firstly, 
excited states undergo a temporal evolution, so that a 
dynamic approach to the CP force may be more appro- 
priate than a static one. Secondly, if an excited atom 
strongly interacts with a body-assisted electromagnetic 
narrow-band field, perturbation theory may break down. 

In the following, we focus on the CP force that results 
from resonant single-quantum exchange between an atom 
and a body-assisted field, which may be described by ap- 
proximating the atom by a two-level system and treat- 
ing the atom-field interaction in rotating-wave approx- 
imation. We will first (Sec. IIII A| develop a simplified 
approximate theory that is similar in spirit with both 
Casimir and Polder's original work Q and the quasi- 
stationary dressed-state approach to strong atom-field 
coupling [fl ||| . We will then (Sec. IIII B)) develop a more 
complete dynamic theory, by starting from the operator- 
valued Lorentz force [lj, [t5| and evaluating its time- 
dependent expectation value by solving the Schrodinger 
equation for chosen initial preparation, thereby taking 
into account the finite bandwidth of a realistic body- 
assisted electromagnetic field. 

A. Static approximation 

As already mentioned at the end of Sec. IIIB1 the 
Hilbert space for the case where a two-level atom that 
is initially prepared in the upper state resonantly in- 
teracts with a body-assisted electromagnetic field can 
be spanned by the states |1)|{0}) and |0)|r A ,w). Using 
Eq. dU [together with Eqs. ||Sj), HD, and @] and re- 
calling the definition (|3"2"|) and the commutation relations 
(127)) and ([55]). we derive 

H\l)\{0}) = E 1 \l)\{0}) 

+ h dwff(r A) w)|0)|rA,w), (39) 
Jo 

fl"|0)|r A ,w) = {Eo + nu)\0)\T At u) 

+ hg(r A ,uj)\l)\{0}). (40) 

Let us assume that the body-assisted field that reso- 
nantly interacts with the atom can be approximated 
by a single Lorentzian-type nonmonochromatic mode v 
of mid-frequency lo v and full width at half maximum 
(FWHM) lv , 
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such that 7„ — > and g 2 (r,u u ) — > oo, but the product 
7 1/ g 2 (r, w„) remains finite, thus 



(42) 



5 2 (r,w)-> |7T7^3 2 (r,cJ y )5(w - u;„) 

in the limiting case of a monochromatic mode. 

The single-resonance assumption implies that Eq. 
takes the form 



iT|l>|{0}) = £i|l)|{0}) + ±M7 R (r A )|0)|l„), (43) 
where 

n R {r A ) = v /2TT- fv g 2 (r A ,uj u ) (44) 

is the vacuum Rabi frequency, and the state |1„) is de- 
fined by 



tin; 



-Sui/2 x/(u-uj u ) 2 +7 2 /4 



kA,w>, 



where <5w is a measure of the distance between two neigh- 
boring lines (7„<§;(W). On recalling Eq. (pM)) , it is not dif- 
ficult to see that (for 7„ — * 0) it is normalized to unity, 



11M = !• 



(46) 



Within the approximation scheme used, from Eq. (|40p it 
follows that 

ff|0)|l„) = (E Q + n^)\0)\l v ) + |Ml R (r A )|l)|{0}). (47) 

Hence we are left with Eqs. (|43|) and (|47[) revealing that 
the subspace spanned by the states |1)|{0}) and |0)|l, y ) 
is invariant under H which, on this subspace, takes the 
well-known Jaynes-Cummings form 



H, 



a/3 



( e 1 \m K {v A ) \ 

( ifin R (r A ) E +hu v J 



(48) 



(for the level scheme, see Fig. [T]). 

Obviously, the single-resonance assumption together 
with the approximation 7„ — * implies that the cou- 
pling of the atom to the residual continuum of the body- 
assisted field is fully ignored, i.e., both radiative and non- 
radiative losses are disregarded. This means that the re- 
sults found in this way are only valid on time scales which 
are short with respect to the associated relaxation times, 
as will explicitly be shown in Sec. IIII Bl 

Straightforward diagonalization of the Jaynes-Cum- 
mings Hamiltonian (|4"8"]l yields the two eigenenergies 



E± = UEo +E l + huo v ) ± |ttl(r A ) 



with 



where 



fi(r A ) = Jfi R (r A ) + A 2 



(49) 



(50) 



(51) 




E x 



frhJ v 



E 



FIG. 1: Level scheme of a two-level atom interacting with a 
single-resonance electromagnetic field. 



is the atom-field detuning, and the eigenstates (the 
dressed states) can be written in the form 

|+) =cos[0 c (rA)]|l)|{O})+sin[0 c (r A )]|O)|l v ), (52) 
|-) = -sin[0 c (r A )]|l)|{O})+cos[0 c (r A )]|O)|l,), (53) 

where, according to 

tan[20 c (r A )] = -^M ; 9 c (r A ) G [0, 7r/2], (54) 

the coupling angle c (r A ) has been introduced, and the 
relation 



1 



yjl + cot 2 (a) 



= sin(a), a <G [0, 7r] (55) 



has been used. Note that in the case of exact resonance, 
A = 0, Eq. |J5DJ reduces to 

fi(r A ) = n R (r A ), (56) 

and Eqs. ([B^|) and (f53"| simplify to 

|±) = -J=(±|1)|{0}) + 10)11,)). (57) 

Comparing Eq. (|4"9")l with Eq. (|3T|) . we can conclude 
that for a system prepared in the state |+) or |— ), re- 
spectively, the CP potential f/+(r A ) or C/_(r A ) is given 
by 

0±(r A ) = ±|ftfi(r A ), (58) 

from which the associated CP force 

F±(r A ) = -V A C/±(r A ) (59) 

can be obtained. Equation generalizes the result de- 
rived in Ref. [5[ for the case of a two-level atom strongly 
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coupled to a single standing wave in a cavity to an arbi- 
trary resonator-like arrangement of bodies giving rise to 
strong atom-field coupling. As can be seen from Eq. (|50|) 
together with Eqs. and (J25J), the CP potential can 

be quite generally expressed in terms of the imaginary 
part of the Green tensor of the bodies as 

U±(r A ) = 

± i V 2 Vo7^^ d io-IniG(r A ,r A ,^)-dio + h 2 A 2 . 

(60) 

At this point we recall that the ground-state variance 
of the electric field in a frequency interval Suj is given 
by [15| 



<[AE(r)] 2 ) fc = [({0}|E 2 (r)|{0}) - ({0}|E(r)|{0}) 2 

= ^t° [ dc^ 2 Im[TrG(r,r,w)], (61) 

which in the single-mode case considered here reduces to 

([AE(r)] 2 ) fc = iVo7^Im[TrG(r,r,^)]. (62) 

Comparing Eq. (J62J) with Eq. ([58]) [recall Eqs. (J25j) and 
(|44[)]. we see that the CP potential is essentially deter- 
mined by the ground-state fluctuations of the electric 
field at the position of the atom along the direction of 
its transition dipole moment. Moreover, we see that for 
a system prepared in the state |+) (|— )), the atom is 
repelled from (attracted towards) regions of high field 
fluctuations. This is in agreement with the result found 
for the case of a standing wave in a cavity [5j , where the 
atom is repelled from (attracted towards) the antinodes 
when the system is in the state |+) (|— )). 

To make contact with the results of perturbation the- 
ory, which are applicable in the case of weak atom-field 
coupling, let us consider the limiting case of large de- 
tuning, |A| ^>r2rt(rA), where, according to Eq. ([54"]) . the 
coupling angle approaches 8 c (r A ) ~ tt/2 and 9 c (v A ) = 0, 
respectively, for positive and negative detuning, and the 
dressed states (152"]) and (f5"3"|) approximate to 



|0)|1„) forA>0, 

-|1)|{0}) forA<0, 

|1)|{0}) forA>0, 

|0)|1„> forA<0. 



(63) 
(64) 



By expanding the square root in Eq. (|50[) and recalling 
Eq. ([37]) . one finds that the CP potentials Ui[ -j(r A ) and 
£^oi„( r A) associated with the states |1)|{0}) and |0)|lj,), 
respectively, can be written as 



^i{o}(r A ) = -—n%_(r A ), 



f/ ffl „(r A ),i^(r A ). 



(65) 
(66) 



Next, we employ the Kramers-Kronig relations for the 
(spectral) response function uj 2 G^(y, r, uj) to note that, 
according to the single-resonance approximation (|42|) 
made, we may write [recall Eqs. (|25p and (|4"4"|) ] 

^ w 2 di • Re G (1) (r, r, uj) ■ d i 
do;' 



,'2 



, to' — uj 
\irh'y l ,g 2 (v A ,u} v )'P 



di -ImG (1) (r, r,u')-d 01 
— d{uj - us v ) 

-co W ~ w 

= 77 RA \ foruj^ujv 67 
i(uj u -u) 

(V: principal value), which is valid up to a position- 
independent constant. Here, the Green tensor has been 
decomposed into the translationally invariant bulk part 
G^ -* and the scattering part 

G(r, r', w) = G (0) (r, r', w) + G (1) (r, r', w). (68) 

Upon discarding position-independent terms, Eqs. (|65p 
and (|66[) can hence be rewritten in the form 

[7x{ }(r A ) = -Mo^iodio-ReG (1) (rA,rA, ^io)-d i, (69) 
Uoi„(r A ) = Mow 2 odio-ReG (1) (rA,r A , wio)-d i. (70) 

As expected, J7i{o}(i"a), Eq. (|69|) . is just the resonant 
part of the CP potential obtained in leading-order per- 
turbation theory for a two-level atom that is prepared in 
the upper state and (weakly) coupled to a body-assisted 
electromagnetic ground state, cf. Refs. [Hill. The fact 
that only the resonant part of the potential appears is 
obviously due to the rotating-wave approximation. The 
result ([70]) is new, Uoi, y (r A ) gives the resonant part of 
the (weak-coupling) CP potential for the case where in- 
stead of the atom being excited and the field being in 
its ground state, the atom is in its ground state, but the 
field is in single-photon Fock state ([45]) . It is seen that 
the potentials Uoi u (r A ) and Ui{ y(r A ) and hence also the 
associated forces just differ in sign. Regarding the forces 
as being the result of recoil due to emission/ absorption of 
a photon (where in the presence of bodies the photon is 
predominantly emitted/absorbed in a certain direction, 
leading to a nonvanishing net recoil) , this difference obvi- 
ously reflects the fact that for the excited atom the reso- 
nant process responsible for the potential is the emission 
of a photon while for the field being excited the resonant 
process is the absorption of a photon (hence the direc- 
tion of net recoil is reversed with respect to the emis- 
sion case). It should be pointed out that the validity of 
Eq. (|69p does not require a single-mode approximation, 
so this equation is generally valid in the weak-coupling 
regime. Clearly, the state |1)|{0}) can then no longer be 
regarded as being a (quasi-)stationary one in general. 

Finally, let us consider the situation where the system 
is prepared in a superposition state 



= cos0|l)|{O}> 



(71) 
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9 e [0, 7r], which includes the special cases 



and 



|0=O) = |1)|{O}), 
|0=tt/2) = |0)|1„ 

|0 = 0c> = |+>, 

\9 = c + ir/2) = |- 



(72) 
(73) 
(74) 
(75) 



Note that when the system is prepared in a state \6) 
whose projections onto both |+) and |— ) do not vanish, 
Rabi oscillations will be observed, leading to a nontrivial 
dynamics of the CP potential and the associated force, as 
we will show in Sec. 111113] Here, we restrict our attention 
to the potential and the force at initial time, where the 
system is in the state \9). 

The CP potential as the position-dependent part of the 
expectation value of the Jaynes-Cummings Hamiltonian 
of the system prepared in a state \9) is calculated to be 

Ue(r A ) = \(9\+)fU + (r A ) + \(9\-)fU-(r A ) 
= cos 2 [9-9 c (r A )]U+(r A ) + sin 2 [9-9 c (r A )}U^(r A ) 

= ±hcos{2[9-9 c {r A )]}n(r A ), (76) 



where Eqs. (j52|) , ([53]) and (|5S|) have been used. According 
to Eq. (l38|) . the associated CP force then follows as 

F e (r A ) = - ±/icos{2[0-0 c (r A )]}V A fi(r A ) 

- hsm{2[9~9 c {r A )]}Q(r A )V A 9 c (r A ). (77) 



In order to carry out the derivatives in Eq. (I77|) . we note 
that Eqs. §5U§, (O and imply 

sin[20 c (r A )]fi(r A ) - n R (r A ), (78) 
V A 0(r A ) = sin[20 c (r A )]V A fi R (r A ), (79) 

while taking the derivative of Eq. (|54|) and using the re- 
lations 



1 J cos(a) for a € [0, 7r/2], 

+ tan 2 (a) [-cos(a) for a G [tt/2, tt], 

we derive 



V A c (r A ) 



cos 2 [2fl e (r A )]V A f! R (r A ) 
2A 



(80) 



(81) 



Using Eqs. (J5J), (J7SJ), ([75]). and ([HI]), we can eventually 
write Eq. (|TT[) in the form 



F e (r A ) = -\h (cos{2[0-0 c (r A )]} + cot[2# c (r A )] 

xsin{2[0-0 c (r A )]}) V A 0(r A ), (82) 

which of course contains the previous result (|59p in the 
special case of the system being prepared in one of the 
dressed states |±), recall Eqs. (J74j) and (|75)) . 

In the case of exact resonance, A = 0, where 9 C — tt/4, 
Eqs. ([75]) and ^ reduce to 



F e (r A ) = -iftsm(20)V A Q R (r A ), (84) 

respectively. In complete analogy to the derivation of 
Eqs. (f6"9"|) and ([7D|) . one can show that for large detun- 
ing, |A|;^r2 R (r A ), which corresponds to weak atom-field 
coupling, Eqs. ([75]) and (|52"j) approximate to 



r„lr.v) : -A cos(20)tt R (r A ) 



cos(26»)^ w 2 di -Re G (1) (r A , r A , wi )-d x (85) 



and 



F e (r A ) = — cos(20)ft R (r A )V A f> R (r A ) 

= cos(20)/i oW 2 o V A d lo -ReG( 1 Hr A ,r A ,w lo )-d ol , (86) 

respectively. 

B. Dynamical theory 

In order to develop a dynamical theory of the CP force, 
which particularly takes into account the line width of the 
nonmonochromatic mode interacting with the atom, we 
start from the operator-valued Lorentz force acting on 
the atom in electric-dipole approximation [14j . 



F=<| V[d-E(r)]+-[dxB(r)] 



pel _|_ pmag 



(87) 



where we have separated the force into its electric and 
magnetic parts. Note that while we have expressed the 
force in terms of E(r) which has the physical meaning 
of a displacement field, Eq. (|87j) is also valid with the 
(physical) electric field in place of E(r) [3, [3]. For a two- 
level atom whose interaction with the electromagnetic 
field is treated within the rotating-wave approximation, 
we may employ Eqs. (JTTJl, (22), and to present F el 
in the form 



-h I duj 
'a 



Vg{v,u>)a(Y,u)a ] 



H.c. 



(88) 



and by using Eq. (fl"4"|) , we may write F mag as 

k A— e,m 

P°° A 1 

x/ —d 10 x[VxG x (r,r',u;)}f x (r',^} +H.c. 
Jo w \ . 

' r— r a 



U e (r A ) = ifisin(20)r! R (r A ) 



(83) 



(89) 
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For chosen atomic position (r A *—>■ ta), the dynamical 
CP force is just the average Lorentz force 

F(r A ,t) = (m\F cl \m) + (m\F mag \m) 

= F c \r A ,t)+F^(r A ,t), (90) 

where the state of the system, \ip(t)), evolves according 
to the Schrodinger equation 



ih2 j; W)) = H\m)> 



(91) 



with the Hamiltonian being given by Eq. |T]). together 
with Eqs. ©, dHJ, and @. 



1. Quantum- state evolution 

In order to include in the evolution of the quantum 
state the loss mechanisms associated with the line width 
of the nonmonochromatic mode that strongly interacts 
with the atom as well as the interaction of the atom 
with the residual, weakly interacting field, we approxi- 
mate g 2 (r, uj), Eq. ([25]) , as 



5 2 (r,w) = gl(r,u v ) 



(w - UJ y ) 2 + 7 2/4 



+ <?' 2 (r,w), (92) 



with the first and second term corresponding to the con- 
tributions from mode v and the residual field, respec- 
tively. In contrast to the drastic approximation that we 
have made in Sec. IIII Al [recall Eq. (|4"2f ], the width 7„ of 
the mode is now finite, but it should still be small com- 
pared to the separation of the mode from the neighboring 
ones, 7^ <C Su>. Furthermore, we have retained the resid- 
ual field continuum g' 2 (r, uj) , which is assumed to be a 
slowly varying function of uj in the vicinity of uj = uj v such 
that its effects on the dynamics is adequately described 
by the Markov approximation, and 



u u +Suj/2 
uj v — 8uj/2 



duj g' 2 (r,u) < / dtu 



„-Su>/2 {u-U v ) 2 +^l/A 

= {n 2 R {r A ) forr = r A . (93) 

Under these assumptions, the natural generalization of 
the state \l u ) defined in Eq. (|45p reads 



1 



tt 7i/ g{r A ,u;„) J Uu _ 5u}/ 



u„+Suj/2 



&u>g(r A ,uj)\Y A , uj). 



(94) 

This state is normalized to unity, as can be easily seen 
by using Eqs. ((34|) and (l93|) : it reduces to the state ([45]) 
when setting 7„ — > and g' 2 (r, uj) — > 0. 

Let the system at initial time to be prepared in a state 
\ip(t )} = \6) , where \6) is defined according to Eq. (fTT) . 
but now with ll^) as given by Eq. (|93|) . Within the 
rotating-wave approximation, the state \ipg(t)) can then 
be expanded as 

^(*)>=V>i(t)|l)|{0})+ / du;V>oM)|0>|rA,u/>, (95) 
Jo 



where the coefficients ipi(t) and i(>o(uj, t) satisfy the initial 
conditions 



^i(to) = cos 9, 



(96) 



M^M) = \f^- 9< f ,U) \ Bind[0(w-w v + *w/2) 
V t^Iu g{r,w v ) 

— <3(uj — LO v — Suj/2)] (97) 

[0(x): unit step function]. We insert Eq. (195)) in Eq. (|9T|) . 
make use of the definition 1|32|). and recall the commu- 
tation relations (|27|) and (|28|) to obtain the system of 
coupled equations 

POO 

i, 1 (t) = -ih~ 1 E 1 ^ 1 {t)-i dujg(r A ,u;)Mu,t), (98) 



^o(w, t) = -i(h 1 E + lu)M^, t) - ig(r A , w)^i(t), (99) 

where ipo(uj,t) can be eliminated by formally solving 
Eq. (|99|) [together with the initial condition (|97|) ]. leading 
to 



V ?nv g(r,uj u ) 
[8(w -w„ + <W/2) - 8(cj - w„ - <W/2)] 

-ig(r A ,u) f dre-^^+^-^Mr)- (100) 

After substitution into Eq. (|9"8")) . the frequency integral 
resulting from the first term in Eq. (|100p can be carried 
out according to 



2 i sin ( 



dw 5 2 (r A ,a;)e- i (- Eo / R + a ')(*- to ) 



7T7„ g(r,u u )J u „_ guj/2 
-i l ^ R (r A )sin0 e [- J ^/ ri+ ^)-^/ 2 ]( t -* ) 



^j(r A ) 
7rfi R (r A ) 



-i(So/»+««)(t-tq) 5l/au ^ _ t() ) ( 101 ) 



[^i/5w(i ~ ^o) : coarse-grained delta function (time scale 
1/Sw)], where we have used Eq. (|9"2")) together with the 
(approximately valid) relation ( 7 „ <C &<;„) 



2?r 



2 ^ (W ~ W„) 2 + 7 2/4 7 „ 



(102) 

The second term in Eq. ([TUT]) , in which r£ (r A ), Eq. ([T08|) . 
is the spontaneous decay rate of the upper atomic state 
due to the weak coupling of the atom to the residual 
background field, is only non- vanishing near t — to, where 
its contribution is negligible in comparison to that of the 
first term [recall Eq. ([9"3"]l]. Hence Eqs. J5S]) and (TODD 
lead to the following closed equation for (i) : 

- |iQ R (r A )sin0 e [-i(*Vft+",)-7,/2](t-to) 



dw.g 2 (r A ,w) / dre' l(£o/fi+w)(t - T) ViM- (103) 

11 J t 
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In order to solve this integro-differential equation, we 
make use of Eq. (f9"2"| and treat the term g' 2 (r,uj), which 
describes the effect of the residual field continuum, on 
the basis of the Markov approximation. It can then be 
shown that (see App. [XJ 

fa(t) = e i-iEl/h-t8u J ' 1 (r A )-r> 1 (r A )/2](t-t a ) ( j )i ( t ^ 

where <j>\ (t) is the solution to the differential equation 

Mt) + {iA(r A ) + [ lu - T[(r A )}/2}Mt) 

+ iO|(r A )^(i) = (105) 

together with the initial conditions 



2. Average Lorentz force 



ii(io) = cos0, 0i (to) = --^ R (r A )sin0. (106) 



Here, 



^o^, / ~j, ., .2 dio-ImG (1) (r A ,r A ,a;)-doi 



SuA(T A ) = / dujLU 



wio(rA) - tJ 



fi*(r A ) A(r A ) 



4 A2(r A )+ 7 , 2 /4 



(107) 



and 



ri(r A ) 



^[wio(r A )] 2 d i ■ ImG[r A ,r A , wi (r A )] ■ d 
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^R^a) Ju 



4 A2(r A )+ 7 2/4 



, (108) 



respectively, are the shift and width of the upper level 
associated with the residual field, and 

wio(r A ) - h- 1 [£i(r A ) - Eq] = uj 10 + 6u[(t a ) (109) 

and 

A(r A ) =w„-Cio(r A ) = A - <L^(r A ) (HO) 

are the respective shifted atomic transition frequency and 
detuning. Note that in Eq. (|107[) . the vacuum Lamb shift 
contribution to the level shift has been absorbed in the 
bare transition frequency uio by making the replacement 
G i ► (?W . Writing the general solution to the differential 
equation (|105[) in the form 

Mt) = c + (r A )e n + (rA)(t ~ to) +c_(r A )e n -( rA ) ( *- t0 \ 

(111) 

we find that 

Ofc(r A ) = -H iA ( r A) + [7,-r'i(r A )]/2} 



T i V /{ lA ( r A) + [7,"r' 1 (r A )]/2} z - Ql(r A ) , (112 
and the initial conditions (I106|) imply 

n=p(r A ) cos6> + ^i£ln(r A ) sm6 



c±(r A ) 



n T (r A ) - fi±(r A ) 



113) 



Substituting Eqs. (88]), (EH), and J95J) into Eq. ([90]), re- 
calling Eq. (|32|) . and using the commutation relations |6]), 
([7]), lj2"7j). and l(2"8j) . we may express the average Lorentz 
force in terms of -0i (t) and tj)o(u>,t) to obtain 



F el (r A ,i) = F c e 1 (r A ,t) = 

r°° [V5 (1) (r,r A ,w) 

-ft/ da>^(t)Vto(w,*)- ; r 

Jo 5(rA,w) 



.c. 

(114) 



and 

F^(r A ,()EF7(r A ,f) 

ITT J dt 

dmx [VxImGW(r,r A ,w)-d i] r= 



9(rA,w) 



+ C.c, (115) 



where g^(r,r' ,uj) is defined according to Eq. (|29[) with 
G(r, r',w) being replaced with G^ 1 - 1 (r, r', uj). Note that 
the bulk part of the Green tensor, (r, r',uj), does not 
contribute to the force. By making use of Eq. (|100[) . we 
can express P| (r A ,t) and F™ ag (r A ,t) entirely in terms 
of ipi(t). For example, Ff(r A , t) can be given in the form 

Ff (r A ,t) = _ ^g e [-W^)-7./a](t-*.)^( f) 
SiR(r A ) 



/>00 

: [Vff« (r, r A , + ih / du [V ff « (r, r A , 

Jo 



OJ 



x / dTe- i ^/* +w ^ t - T Vi(*)^i(T) +C.c, (116) 

J t 

where, in close analogy to the derivation of Eq. Q103p , the 
first term has been obtained by means of the relation 

7,74 



g{r,r',uj) = g v (r,r',uj u ) 



(u - u v f + 7 2/4 



(117) 

[which reduces to Eq. ([92]) for r=r'], upon using Eq. (|102[) 
and discarding the term proportional to 5i/s^{t — to), 
recall Eq. (flQl)) . 

Now we can make use of ipi(t) as given by Eq. (|104|) 
together with Eq. (Tl~TTj) . From Eq. (fTT6|) for Ff(r A ,t) 
and from the analogous equation for F™ ag (r a , t) we then 
obtain, after some algebra, the following expressions for 
the electric and magnetic components of the resonant 
part of the force observed in the case when the system is 
initially prepared in a state of the type ([71]) : 

+ h duo [V<7 (1) (r,rA,w)] r=rA s(r A ,cj,t-io) + C.c., 



(118) 
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F e S (rA,*) = -^y S1 n^-g(r A) t-* ) 
x di x [VxImG^^rAjW^-do^ 



dujLud w x [VxImG^fr, r A ,u;)-doi] 



where 



q(r A ,t) = e 



= {-iA(r A )-[ 7 „+ri(r A )]/2}t 



X — «(r A ,w,t-to) + C.c., (119) 

di and 



C ;(r A )e n + (rA)t + c*_(r A )e n -^ A)t ] (120) 



s(r A ,w,t) =|c+(r A ) 



+ c+(r A )c_(r A ) 
+ cl(r A )c + (r A ) 
+ k-(r A )| 2 



e [-r' 1 (r A )+0;(r A )+^ + (r A )]t _ e {ip' w (r A ) -w] -V\ (r A )/2+n* (r A )}t 

w - £ 10 (r A ) + iri(r A )/2 - ifi+(r A ) 
e [-ri(r A )+n;(r A )+n_(r A )]t e {i[aJ (r A )-w]-ri(r A )/2+n^(r A )}t 



w - c^ (r A ) + »Ti(r A )/2 - iO_(r A ) 

e [-ri(r A )+fT (r A )+fi+(r A )]t _ e {i[£' 10 (r A )-w]-ri(r A )/2+ni(r A )}t 

w - ^ (r A ) + iri(r A )/2 - ifi+(r A ) 
e [-ri( rA )+n* (r A )+n-( rA )]t _ e {»[aio( r A)-w]-ri(r A )/2+«* ( rA )}t 



- w£ (r A ) + iri(r A )/2 - ifi_(r A ) 



r 



(121) 



3. Weak and strong coupling limits 

In order to gain deeper insight into the results, let us 
consider the limiting cases of weak and strong atom-field 
coupling in more detail. To make contact with earlier 
results, we begin with the weak-coupling regime. 

a. Weak atom-field coupling. Weak atom-field cou- 
pling is realized if the width 7„ of the nonmonochromatic 
mode is sufficiently large, 

lv > 2f7 R (r A ), (122) 

and/or the detuning A(r A ) is sufficiently large, 

|A(r A )| » 2f^(r A )/ 7 „. (123) 

In both cases, the first term under the square root in 
Eq. (|112[) is much larger than the second one and a Taylor 
expansion yields the approximations 

n±(r A ) = 

'-»A(r A )- [7,-r;(r A )]/2 
»»|(r A ) A(r A ) 



A 2 (r A )+ 7 2/4 



A2(r A )+ 7 2/4 



(124) 



For a system initially prepared in the state |1)|{0}), i.e., 
6 = [recall Eq. ([72])]. Eq. (fTT5t) approximates to 



c+(r A ) = 0, c_(r A ) = 1 



(125) 



because |0 + (r A )| > |fi_(r A )|. Substituting Eqs. ([103]) . 
(TlTTj) . dim , and (TT23)) into Eq. and recalling 



Eqs. (|107p and (|108[) . we find that in the weak-coupling 
limit the temporal evolution of the system initially pre- 
pared in the state |1)|{0}) may be given by 

\lp(t)) = e [- ,;£ i/ ?i - i ' 5 "i( r A)-ri(r A )/2+0_(r A )](t-to) |1)|{0}) 
— e [-iEl(r A )/ft-ri(r A )/2](t-t )H\| |QT.\ (126) 

Here, the shift and width of the upper atomic level are 
given by 



<5wi(r A ) 



and 



,dio-ImG (1) (r A ,r A ,w)-d i 



wio(r A ) - to 



(127) 



Ti(r A ) = ^[wio(r A )] 2 di TmG[r A ,r A ,cZ'io(r A )] -d i, 

(128) 

respectively, and the shifted atomic transition frequency 
a)io(r A ) reads 

Cj w (r A ) = h- 1 [E 1 {r A )-E ] = uj 10 + <5cji(r A ) (129) 



in place of Eq. Q109p . Note that in the weak-coupling limit 
the two terms in the decomposition (j9"2"|) can be treated 
on an equal footing, so that the decomposition becomes 
superfluous. 

The associated force in the weak-coupling limit can 
most conveniently be derived by returning to Eq. (I116P 
[and the analogous equation for F™ ag (r A , t)] and making 
therein use of Eq. (|126p in the form 

Mt) = e H&(rA)/A-ri(r A )/2]Ct-T) (130) 
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Evaluating the time integral in the spirit of the Markov 
approximation by setting [^(r)! 2 i— ► |^>i(£)| 2 and letting 
the lower integration limit tend to — oo, one finds, on 
recalling Eq. {S3]), that F™ { Q g } (r A , t) = and hence, 

F 1{0} (r A ,<) = Ff {0} (r A) t) = e - r ^)(*-to) Fl{0}(rA)) 

(131) 

where 



Fi {0 }(i-a) 



Mo 



duj w 



[Vdio-ImGWCr.rA.wJ-doi], 



C.c. 



w - u>i (r A ) - ir x (r A )/2 
/z n? (r A ){Vd 10 ■ G«[r,r A ,ft 10 (r A )] ■ d i} r=rA 

+ C.c. (132) 



with 



fiio(r A ) = wi (r A ) + tTi(r A )/2. 



(133) 



As expected, Eq. (|13ip agrees with the resonant contri- 
bution to the force as derived in Refs. [3, [l5|, in the 
special case of a two level atom that is initially prepared 
in the upper state. 

b. Strong atom-field coupling. Having thus estab- 
lished that the general expression (|118|) — (|121|) reproduces 
earlier results in the weak-coupling limit, we now turn our 
attention to the strong-coupling regime. Strong atom- 
field coupling is realized if both the spectral width j u of 
mode v and the width r'^r^) of the upper atomic level, 
which is associated with the residual field, are sufficiently 
small, 



7„,ri(r A ) < 2J2 a (r A ) 



(134) 



and its mid-frequency u) v is sufficiently close to the 
atomic transition frequency cI>io(r A ), 



|A(r A )| « 2fi R (r A )/ 7!/ . 



(135) 



Even when the coupling is at least moderately strong 
such that 7„,r' 1 (r A ) < 2f2;a(r A ) and the inequality (|135|) 
holds [note that in addition, |A(r A )|r£(r A ) < 20|(r A ) 
is automatically fulfilled as a direct consequence of 
Eqs. and (|108[) ], then the real part of the square 
root in Eq. (| 1 1 2[) becomes negligibly small so that 

fi±(r A ) - -H lA ( r A)+[7,-ri(r A )]/2} T ^^(r A ) (136) 

is approximately valid, where f2(r A ) is now given by 



n(r A ) = \/nl(r A ) + A 2 (r A ) - [^-T'^va)] 2 /A (137) 



in place of Eq. ([HO]), Accordingly, ipi(t) [Eq. (fT04| to- 
gether with Eq. (jllip ] approximately reads 

Mt) = e -7(rA)(t- t o)/2L +(rA)e -^ + (r A )(t- t „)M 

+ c_(r A )e- jiJ - (rA)( *- to)/?i l, (138) 



^±£n(rx) (139) 



where 

ii±(r A ) = 



are the eigenenergies of the system in place of Eq. (|49|) . 

7 (r A ) = |[ 7 ,+r' 1 (r A )] (140) 

is the total damping rate and the coefficients c±(r A ) are 
given by Eq. (TT5J) with fi±(r A ) from Eq. (fT36"j) . Note that 
from Eq. ifliD]) [together with Eq. (fTOg]) ] and Eq. (fT07|) 
it follows that 



y v /2, |A(r A )| « 7 „/2, 

ri(r A )/2, j„/2 « |A(r A )| « 2f> 2 (r A )/ 7 „ 

(141) 



fi 2 l (r A )A(r A )/ 7 2 , |A(r A )| « 7 „/2, 

fei(r A ), 7k /2 < |A(r A )| < 2f7 2 ,(r A )/ 7 „ 

(142) 




[recall Eqs. (fl2?|) and (fl28])]. 

We substitute Eq. <p~36|) into Eqs. (fUg]) and (fTT9|) and 
calculate the frequency integrals by means of Eq. (|117p . 
The residual-field term g'(r,r',uj) can be pulled out of 
the integration and does hence not contribute to the res- 
onant force [since the only remaining poles giving rise 
to resonant terms are then those of s(r A ,u>,t — to), at 
which the respective numerators vanish, cf. Eq. (| 1 2 . 
The nonvanishing contribution due to mode v can be 
evaluated in the spirit of Eq. (|102|) . and after a tedious, 
but straightforward calculation, we arrive at 



F e e '(r A ,i) = -h7Tj u [Vg^(r,r A , l j v )] i 



- 7 (r A )(t-t ) 



|c + (r A )| 2 + cl(r A ) C+ (r A )e-' n ^)( t - t °) 
A(r A ) - fi(r A ) -i[7„-ri(r A )]/2 

|c^(r A )| 2 + c;(r A ) C _(r A )e^ r A)( t ^o) 
A(r A ) + fi(r A ) - i[7^-ri(r A )]/2 



+ C.c. 

(143) 



and 



F 6T g ( r A,0 = M0^7^^R( r A) 

x di x [VxImG (1) (r,rA,w„)-d i] 
c*_(r A )c+(r A )e 



- 7 (r A )(*-t ) 



r=r A 
iO(r A )(t-t ) 



A(r A ) — 51(r A ) — i [ lu - T' x (r A )] /2 
C ;(r A )c_(r A )e i0 ( rA )( t -*°) 



A(r A ) + fi(r A ) -i[ju~ n(r A )] /2 



+ C.c. (144) 



In order to compare with the (purely electric) weak- 
coupling result (I131|) - (|133j) . let us consider the case where 
the system is initially prepared in the state |1)|{0}) [i.e., 
6 = 0, recall Eq. {72])] and examine the electric part of the 



12 



force in more detail. Noting that the coefficients c±(r A ), 
Eq. ifTTg]) , reduce to 

C±(rA) = 2^ ' (W5) 

so that for real dipole matrix elements where 



V 3 w (r,r A ,w„)| r=rA = V ff w *(r,r A ,^)|r=r A 

= g(r A ,u v )^ A g(r A ,u v ) (146) 



[recall Eqs. HH}, (JUJ), and (|29|)]. from Eq. 1(135]) it follows 
that 



FiWO-A,*) 



_ p -7(r A )(*-*o) 



sin 2 [r!(r A )(t~to)/2] 



?i7r 7l , [V#W (r, r A , uv)] r=r; 

X A(r A ) + z[ 7t ,-ri(r A )]/2 

A 2 (r A ) - [T^TUrAjjVj 
A2(r A )-[ 7l/ -r' 1 (r A )P/4+0|(r A ) 



+ C.c. (147) 



Since, within the approximation scheme used [Eq. (|117[) ]. 

we may set 

— ^d 10 -G UJ (r,r',w)-d i = — — -, (148) 

fi,7T CJ^ — 0J + l^ v j I 



we may rewrite Eq. (1147)) as 

F? {0} (r A ,t) = 2e-^ rA )(*-*°) sin 2 [fi(r A )(t - t )/2] 



where 
C(r A ) 



xC(r A )F 1{0} (r A ), (149) 
A 2 (r A )-[ 7l/ -ri(r A )/2] 2 /4 



A 2 (r A ) - [ 7 ,-r' 1 (r A )/2]2/4 + fif l (r A ) 
and F 1 { }(r A ) is given by Eq. (I132|) with 

n 10 (r A ) = Q' w (r A ) + *Ti(r A )/2 
in place of Eq. JT55J). Note that 

7^/4 



, (150) 



(151) 



C(r A ) 



7 2 /4-^(r A )' l A ('A)l<7,A 

A 2 (r A ) 
A2(r A )+Q|(r A ) ' 

7 ,/2<|A(r A )|<2f> 2 (r A )/ 7l ,. 



(152) 



Equations (]149|) — (|152[) reveal that the electric part of 
the force for moderately strong to strong atom-field cou- 
pling differs from the respective weak-coupling result 
(|13ip — (|133p in several respects. Firstly, the strength of 
the force is modified by the correction factor C(r A ). 
Secondly, and most strikingly, the time dependence of 
the force is given not by a simple exponential decay, 
but by damped Rabi oscillations with frequency f2(r A ) 



and damping rate 7 (r A ) — oscillations that are well pro- 
nounced when the inequalities (|134[) hold. Note that 
only the magnitude, and not the sign of the force is 
oscillating, with the sign being determined by the sign 
of the detuning A(r A ), as for the force in the case of 
weak atom-field coupling. As seen from Eq. (|14ip . for 
small detuning, |A(r A )| <C 7 „/2, the damping is domi- 
nated by the radiative and non-radiative losses the quasi- 
monochromatic mode suffers from, while radiative and 
non-radiative spontaneous decay of the upper atomic 
state is the dominant loss mechanism for large detuning, 
7 „/2 <C |A(r A )| <C 2f2 2 l (r A )/ 7l/ , where in both limits the 
damping of the force is characterized by one-half the re- 
spective damping rates, 7 „ and Ti(r A ), respectively. This 
can easily be understood from the fact that the force 
depends on the product of the electric-dipole moment 
of the atom and the electric-field strength of the quasi- 
monochromatic mode, whose damping rates are Ti (r A )/2 
and 7„/2, respectively. 

Let us return to the case of arbitrary initial prepara- 
tion of the system. In the limit of strong atom-field cou- 
pling when the inequalities (|134p and (|135|) are fulfilled, 
Eq. (fT3T]) simplifies to 



fi(r A ) = x /r» 2 i (r A ) + A2(r A ) 



l r\ l aj t " a) ■ (153) 
Similarly, after introducing the coupling angle 



tan[20 c (r A )] = 6 c (r A ) G [0, tt/2] (154) 
A(r A ) 

[which replaces Eq. ([54]) ] and using the identities ((53]) and 
(150]) . we may write the coefficients c±(r A ), Eq. (|113p . as 

c+(r A ) =cos 2 [6» c (r A )] cos 6 + sin[0 c (r A )] cos[6» c (r A )] sin0 
= cos[0 c (r A )]cos[0-0 c (r A )], (155) 

c_(r A ) =sin 2 [6> c (r A )] cos6» - sin[6» c (r A )] cos[6» c (r A )] sin6» 
= - sin[0 c (r A )] sin[0 - c (r A )]. (156) 



Equation (|143|) for the electric part of the force then reads 



lr = r A e -7(r A )(t-t ) 



2J2 R (r A ) 

x (sin[2# c (r A )]cos{2[0-0 c (r A )]} + sin{2[#-# c (r A )]} 

+ {cos 2 [0 c (r A )]e^ (rA) ( t -*°)-sin 2 [0 c (r A )] e - lO ( rA )(*- to j) 

+ C.c. (157) 

Here, we have used the identities 
^R.(r A ) 



A(r A ) - fi(r A ) 

^R(r A ) 
A(r A ) + fi(r A ) 



= -tan0 c (r A ), (158) 
= cotfl c (r A ), (159) 
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which follow from Eqs. (| 1 53|) and (|154[) upon applying 
Eq. ([55)) . Using Eq. (| 146[) for real dipole matrix elements, 
recalling Eq. (|44|) and employing the relation 

Vy / Q|W + A2(r A )| r=rA = sin[20 c (r A )]V A r!R(r A ), 

(160) 

as implied by Eqs. ([55]) . (|153p . and (|154p . we arrive at 
the final result 

Ff(r A ,t) = e-^ rA )(*- to )(cos{2[0-^(r A )]} 
+ cot[2# c (r A )] sm{2[9-9 c {r A )}} cos[ft(r A )(t - t )]) 

xF+(r A ), (161) 

where 

F + (r A ) = ~\ V^nl(r) + A^(r A )\ r=rA . (162) 
We first observe that at initial time to, Eq. (|161|) reads 

Ff(r A ,t ) = ~ (cos{2[0-0 c (r A )]} + cot[20 c (r A )] 

x sin{2[0-0 c (r A )]}) V^/fi2 ( r ) + A 2 (r A )| r=rA , (163) 

which agrees with the result ([52"| [together with Eq. ([50)1 ] 
found in Sec. IIII A) when neglecting the frequency shift 
5ui[(r A ), i.e., making the replacement A(r A ) i— > A. Equa- 
tion (|16ip further shows that the electric part of the 
force as a function of time is always damped by an over- 
all exponential factor and that it contains an oscillat- 
ing term whose relative strength depends on both the 
coupling angle 6* c (r A ) and the initial state of the sys- 
tem, 1 9). This is illustrated in Fig. [21 where the time de- 
pendence of the electric part of the force is displayed for 
various initial states \9) and fixed coupling angle 26* c (r A ) 
= 3tt/8. The curves in the figure can be grouped in 
pairs of curves differing only by their sign, each of these 
pairs corresponds to values of 9 which differ by 2A9 = it. 
The existence of such pairs is an obvious consequence of 
Eq. p61[) . The figure further reveals that there are two 
extremes of behavior: While for the initial states with 
29 = 29 c (r A ), 26* c (r A )+7T, the force shows no oscillations 
and is purely exponentially damped as a function of time 
[curves (ai) and (as)], the initial states 29= 29 c (r A )+n/2, 
26> c (r A ) + 37r/2 lead to oscillations of maximal amplitude 
around zero [curves (A2) and (84)]. For other values of 9, 
the temporal behavior of the force is a superposition of 
oscillating and non-oscillating components [curves (bi)- 
(b 4 )]. 

In the case 29 = 0, which corresponds to the initial state 
|1)|{0}), the oscillating and nonoscillating contributions 
to the force combine in such a way that the sign of the 
force remains unchanged for all times [curve (bi)] — in 
agreement with Eq. (|149p . The same is valid for the force 
in the case of the initial state |0)|1„) (29 = ir) which has 
just the opposite global sign [curve (b^)]. At a first glance, 



< 




7T 2f 3f 47T 5lT 67T 




n(r A )(t-t ) 



FIG. 2: The time dependence of the electric part of the res- 
onant CP force for strong atom-field coupling is displayed 
for different values of 6, i.e., 26 = 26> c (r A ) (ai), 20 o (r A )+7r/2 
(a 2 ), 2<9 c (r A )+7r (a 3 ), 20 c (r A )+37r/2 (a 4 ),0 (bi), tt/2 (b 2 ), tt 
(bs), 37r/2 (D4), with parameters 2# c (r A ) = 3tt/8 and 7(r A ) = 
0.05 fl(r A ). The angles 29 for the various curves are indicated 
in the small polar diagram. 



one might expect that the respective states evolve into 
each another during the course of time so that the dif- 
ferent global signs would be contradictory. To resolve the 
apparent contradiction, one must ask whether the initial 
state |1)|{0}) ever evolves into the state |0)|1„). A simple 
calculation reveals that the answer is never: Combining 
Eqs. (fTM)) . (fTTTT) . (fHrJj) . (fT55)) . and (fT56]) for 29 = 0, we 
find that for the state \ip(t)) that is initially prepared in 
the state |1)|{0}) the probability [^(rj)! 2 reads 

\Mt)\ 2 = cos 4 [0 c (r A )] +sin 4 [0 c (r A )] + 2 cos 2 [0 c (r A )] 
x sin 2 [6> c (r A )] cos[ft(r A )(< - *o)] 

> cos 2 [20 c (r A )], (164) 

showing that the state |0)|1„) is indeed never reached 
unless in the case of exact resonance, A(r A ) = [9 c (r A ) 
=7r/4], where the (resonant parts of the) forces associated 
with both initial states identically vanish. 

Using the same approximations leading from Eq. (|143p 
to Eq. I|16ip , we find that in the limit of strong atom-field 
coupling, the magnetic component of the force, Eq. (|144p , 
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can be given as 



(r A ,i) 



-7(r A )(*-to) 



sin{2[#-# c (r A )]} 



cos[Q(r A )(t - t )} 



sin[2# c (r A )] 

x ^ cj l/ 7 t/ d 10 x [VxG^^rAj^-doi] . (165) 



Comparing this equation with Eq. (|161[) for the electric 
component of the force, we see that the magnetic compo- 
nent has a quite different vector structure than the elec- 
tric one, and its order of magnitude is roughly f2(r A )/w„ 
times that of the electric component, so that it might 
become relevant in the context of the recently consid- 
ered superstrong coupling regime [2{|. In particular, the 
magnetic component of the force vanishes when the sys- 
tem is initially prepared in the state \0 = O c (ta)) or the 
state \9 = 6* c (r A ) + 7r/2) for which the electric compo- 
nent of the force is nonoscillating. In all the other cases, 
the magnetic component is — in contrast to the electric 
component — always purely oscillating around zero. 
Finally, let us compare the nonoscillating force 



fi*(r) 



A 2 (r A )| 



(166) 

which is observed when the system is initially prepared 
in an approximate energy eigenstate, i.e., |0 = # G (r A )) or 
\6 = # c (r A ) + 7r/2), with the corresponding static result 
[Eq. pi together with Eqs. ([50|) and ([58])] . We see that 
(i) the static approximation may be regarded as being a 
good approximation on time scales which are small com- 
pared to 7 _1 (r A ), and (ii) the atom-field detuning (|110[) . 
which enters the force, is different from its bare value (|5"Tj) 
due to the coupling with the residual field. However, from 
Eq. (|142|) this effect may be expected to be small in gen- 
eral. 



IV. SUMMARY 

Based on macroscopic QED in linear media, we have 
developed a general theory of the resonant CP force ex- 
perienced by a two-level atom in the presence of arbitrary 
linear bodies, with special emphasis on strong atom-field 
coupling. Assuming that the initial state is a (coher- 
ent) superposition of states that carry a single excitation 
quantum each, we have first worked within a static ap- 
proximation. Reducing the Hilbert space of the system 
to an approximately invariant two-dimensional subspace 
on which the Hamiltonian assumes a Jaynes-Cummings 
form, the eigenenergies and eigenstates have been con- 
structed according to the well-known dressed-states ap- 
proach. Identifying the position-dependent part of the 
eigenenergies with the CP potential for the system being 
prepared in a dressed state, a simple intuitive picture for 
the CP force has been obtained, generalizing results ob- 
tained for cavities to arbitrary resonator-like equipments. 



As the static approximation does not take into account 
the decay of excited states due to unavoidable radia- 
tive and non-radiative losses, the result found for a sys- 
tem prepared in a dressed state is only valid on a time 
scale which is short compared to the time scale of de- 
cay. For systems initially prepared in other than dressed 
states, the static approximation is even more problem- 
atic, because it also neglects the Rabi dynamics which 
is to be expected in the strong-coupling regime. Moti- 
vated by these shortcomings, we have developed an alter- 
native, dynamical approach to the problem by starting 
from the operator-valued Lorentz force and identifying 
the CP force with its expectation value where the respec- 
tive state vector of the system solves the time-dependent 
Schrodinger equation for given initial condition. By sep- 
arating the body-assisted field into the part that is in 
(quasi-)resonance with the atomic transition and strongly 
interacts with the atom and the residual part that weakly 
interacts with the atom, we have solved the Schrodinger 
equation in rotating wave approximation to get a solu- 
tion which fully incorporates the dynamics induced by 
both parts of the field. As a consequence, a general ex- 
pression for the time-dependent resonant CP force has 
been obtained. 



For weak atom-field coupling, this expression repro- 
duces the force obtained earlier for the case where the 
atom is initially in the upper state and the field is in 
the vacuum state. The dynamic behavior of the force, 
which at the initial time effectively agrees with the force 
obtained in leading-order time- independent perturbation 
theory, is given by an exponential damping due to sponta- 
neous decay. For strong atom-field coupling, different dy- 
namical behaviors are possible, depending on the initial 
preparation of the combined system. When the system 
is prepared in a dressed state, the force which initially 
agrees with the force obtained in static approximation, 
undergoes an exponential decay in the further course of 
time which results from the width of the upper atomic 
energy level and the width of the body-assisted non- 
monochromatic mode. When the initial state is a more 
general (single excited) state of the atom-field system, 
then damped Rabi oscillations of the force are observed, 
whose amplitude and mean value sensitively depend on 
the chosen initial state. In particular, when the atom is 
initially excited with the field being in the vacuum state, 
the force due to the electric field exhibits two major dif- 
ferences with respect to the weak-coupling result: It un- 
dergoes damped Rabi oscillations and it is scaled by a 
correction factor. Furthermore, it has been found that 
while the dressed-state force is entirely due to the interac- 
tion of the atom with the body-assisted electric field, for 
a more general initial state, additional oscillating force 
components appear that result from the interaction of 
the atom with the magnetic field. The general results 
obtained can be applied to actual geometries by using 
the appropriate Green tensors, in order to analyze the 
respective spatial structure of the force in more detail. 
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APPENDIX A: SINGLE-RESONANCE 
APPROXIMATION 



Applying Eq. (|92|) . the last term on the r.h.s. of 
Eq. (|103p can be written as 



/■oo />£ 

/ dujg 2 (r A ,uj) / d 

JO J t 



-i(E /h+u)(t-T) 



Mr) 



dLog' 2 {r A ,Lu) / dre- < ( Bo / fi+a, »*- T Vi(r) 

J t 

4 J (cj -u v y + 7 2 /4 

x f dTe~ l( - E " /h+uj)(t - T} iP 1 (T), (Al) 

where, in accordance with the assumptions of the single 
resonance approximation, the first term may be treated 
within the Markov approximation. We hence assume that 
the function e lEl ^ CA -^ t / h, ipi{t) is a slowly varying function, 
so that 



duog' 2 {v A ,uj) / dre-^^+^'^ViM 

I] J t 

I* OO r-t 

Mt) / dLug' 2 (r A ,uj) / dr 
JO J t 



x e 



-i\hcj-E 1 (r A )+E \(t-T)/H 



=i-Vi(t) / d^ ff ' 2 (r A , W )C([^i(r A )- J Bo]A-w) 
Jo 

(A2) 

[C(x) =tt(5(.t) Using Eqs. J35J and ((92) as well 

as the identity 



V 



1 



„ [£?i(r A )-£\,]A-w (^-^) 2 + 7 2/4 
2tt [^i(r A )--Eo]/fi-^ 

{[^ 1 (r A )- J B ]A-^} 2 + 7 2 /4' 



(A3) 



which can easily be verified by means of contour-integral 
techniques for r y u /2^uj u , we find 



dug' 2 (r A ,u) / dre-'^/^^ViW 



-^(rAj-iri^A) ^(t), (A4) 



where we have introduced Eqs. (|10T[> (|110|) . Substituting 
this back into Eq. (|A1|) and evaluating the second term 
by means of Eqs. (|4"4~)) and (|102p . we arrive at 



dug 2 (r A ,Lu) / dre-^^+^-^Mr) 



-«wi(r A ) " |ri(r A )JViW ~ Jfi£(r A ) 
x ^ dTeHW^^M^^), (A5) 

J to 



so Eq. p03|) can be written in the form 



Mt) 



E 1 



-i^L_i^(r A )_ |ri(r A ) 



^i(i) 

^0 R (r A )sin6»e^ l ( Bo /' i +^)-^/ 2 ]( t - to ) 



±Q R (r A ) / dret-^^/^-)-^/ 2 ^*-^^^). (A6) 

J to 



By using Eq. (|104[) , this result is transformed to 



Mt) = -i^ R (r A )sin6»e^ 4A ( rA ^^- r 'i( rA )]/ 2 >( t - t °) 
- ±Q R (r A ) f'dre<- ,A ( rA )-^- r ;( rA )l/ 2 >< i -^ 1 (r) 

■/to 

(A7) 



and after differentiating w.r.t. t, we arrive at Eq. (|105p 
together with Eq. (flTigj) . 
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